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Halos at 22° from the sun attributed to randomly-orientated, pristine hexagonal crystals 
are frequently observed through ice clouds. These frequent sightings of halos formed by 
pristine crystals pose an apparent inconsistency with the dominance of distorted, non- 
pristine ice crystals indicated by in situ and remote sensing data. Furthermore, the 46° 
halo, which is associated with pristine hexagonal crystals as well, is observed far less 
frequently than the 22° halo. Considering that plausible mechanisms that could cause 
crystal distortion such as aggregation, sublimation, riming and collisions are stochastic 
processes that likely lead to distributions of crystals with varying distortion levels, here 
the presence of the 22° and 46° halo features in phase functions of mixtures of pristine 
and distorted hexagonal ice crystals is examined. We conclude that the 22° halo feature is 
generally present if the contribution by pristine crystals to the total scattering cross 
section is greater than only about 10% in the case of compact particles or columns, and 
greater than about 40% for plates. The 46° halo feature is present only if the mean 
distortion level is low and the contribution of pristine crystals to the total scattering cross 
section is above about 20%, 50% and 70%, in the case of compact crystals, plates and 
columns, respectively. These results indicate that frequent sightings of 22° halos are not 
inconsistent with the observed dominance of distorted, non-pristine ice crystals. Further- 
more, the low mean distortion levels and large contributions by pristine crystals needed 
to produce the 46° halo features provide a potential explanation of the common sighting 
of the 22° halo without any detectable 46° halo. 

© 2014 Elsevier Ltd. All rights reserved. 


1. Introduction 

Halos at 22° from the sun are among the most fre- 
quently observed atmospheric optical phenomena [1-3]. 
For example, Sassen et al. [1] reported that 37.3% of the 
daytime sky observations in their ~ 10 years record over 
Salt Lake City, Utah, showed indications of the 22° halo, 
with bright and prolonged halos occurring in 6% of the 
record. Such halos are attributed to randomly-orientated 
hexagonal ice crystals with size parameters of the order of 
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about 100 or greater in cirrus clouds [4,5]. Other less 
frequent optical phenomena, such as the 46° halo, are 
known to be caused by such hexagonal crystals as well. 
It is well understood from geometric optics how refraction 
of light rays through pristine hexagonal ice prisms with 
random orientation forms such halos [e.g., 6,2]. However, 
in situ observations in natural ice clouds indicate a general 
dominance of distorted, aggregated, non-pristine ice crys- 
tals [7-14]. Furthermore, satellite observations that sample 
the total and polarized light reflected towards back- 
scattering angles also indicate non-pristine ice crystals, 
possibly with microscale surface roughness, to be predo- 
minant in tops of ice clouds [e.g., 15-18]. This distortion 
and roughening of ice crystals leads to a randomization of 
the refraction angles between crystal facets, which in turn 
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generally yields featureless phase functions and therefore 
a lack of halos [6,19]. Thus, the frequent sightings of halos 
on the one hand and dominance of distorted, non-pristine 
ice crystals indicated by in situ and remote sensing data on 
the other hand pose an apparent inconsistency. Also the 
huge difference in frequency of sightings of the 22° and 
46° halo [1,3], both formed by randomly oriented pristine 
hexagonal ice crystals, has been interest of past research. 
For instance, Shcherbakov [19] argued that the common 
situation of the sighting of the 22° halo without the 
presence of the 46° halo may be due to the anisotropic 
roughness of the prismatic facets of hexagonal ice crystals. 
The anisotropic roughness model may also be a plausible 
proxy for ensembles of pristine, hollow-ended columns 
and/or bullets [20]. 

Apart from diminishing the halos, distortion or rough- 
ening of hexagonal ice prisms also leads to a significant 
lowering of the asymmetry parameter of their scattering 
phase function [6,21-23]. Thus, the degree to which 
natural cirrus clouds contain pristine particles may have 
important implications for their influence on global cli- 
mate [21,24-26]. 

Although the observability of halos in cirrus clouds may 
depend on cloud parameters such as optical thickness and 
extinction profile, the primary requirement for halos to be 
observable is that the halo features are present in the ice 
crystal scattering phase functions. Usually, studies on the 
presence of halo features in phase functions and how they 
relate to crystal distortion only consider collections of ice 
crystals with a single degree of distortion. However, in situ 
observations generally show a wide variety of crystal 
shapes in a single cloud volume [e.g., 12]. Furthermore, 
possible mechanisms that could cause crystal distortion 
and surface roughness such as aggregation, sublimation, 
riming and collisions are stochastic processes that likely 
lead to distributions of crystals with varying distortion 
levels [9,27]. In this paper, we explore the presence of the 
22° and 46° halo features in phase functions of collections 
of ice crystals with varying distributions of distortion 
levels. 

We describe the simulations in Section 2, show results 
in Section 3 and conclude in Section 4. 

2. Simulations 

Ice crystal phase functions are calculated using the 
standard geometric optics code developed by Macke 
et al. [6]. This ray tracing code takes distortion of ice 
crystals into account in a statistical manner by perturbing, 
for each interaction with a ray, the normal of the crystal 
surface from its nominal orientation by an angle varied 
randomly with uniform distribution between 0° and 
5 x 90°, where 5 is referred to as the distortion parameter. 
Thus, this approach represents the stochastic large-scale 
distortion of a collection of ice crystals. Liu et al. [23] found 
that such an approach is also an efficient, yet relatively 
accurate treatment of low to moderate microscale surface 
roughness, although for more severe roughness some 
significant differences with exact calculations were found. 
For a large collection of ice crystals, microscale surface 
roughness and large-scale particle distortion both lead to a 


similar randomization of the angles between crystal facets, 
which in turn leads to the diminishing of features in the 
scattering phase matrix. Increasing the number of impu- 
rities within ice crystals also has a similar effect [28]. Thus, 
we consider the distortion used here as a proxy of 
randomization of the angles between crystal facets possi- 
bly caused by any of these effects [cf. 27]. Although several 
arguably more realistic, alternatives to the uniform dis- 
tribution of crystal facet tilt angles have been proposed 
[29,30,19,27], the general conclusions of this work are 
considered to be independent of the distortion distribution 
shape, as further argued below. 

Calculations are performed at a wavelength of 555 nm 
with a refractive index of 1.3108 + i2.564 x 10“ 9 [31]. The 
distortion parameter is varied between 0 and 0.8 in steps 
of 0.05 and results are linearly interpolated in between. 
The aspect ratio of columns is varied between 1 and 50 
with 26 geometrically increasing steps. The aspect ratios of 
plates are the inverse of those for columns, for a total of 51 
aspect ratios. Similarly as described by van Diedenhoven 
et al. [22], sizes are varied so that the projected areas of 
the particles, assuming random orientation, correspond to 
the projected areas of spheres with radii of 40, 56, 80, 113, 
160, 226 and 320 pm. This leads to maximum and effective 
crystal dimensions that vary with aspect ratio, but since 
these calculations are performed at a wavelength where 
ice absorption is very weak, the geometric optics calcula- 
tions are size invariant. The phase function calculations are 
performed on a 1 -degree resolution and interpolated in 
between. 

As criteria for the presence of the halo features, we use 
halo ratios as defined by Gayet et al. [32] and Shcherbakov 
[19], i.e., 


where P(6> i) and P(64) are values of the scattering phase 
function evaluated at scattering angles 64 and 0 2 , respec- 
tively. Following Shcherbakov [19], for the 22° halo ratio 
(h 22 ) we use 64=22° and 0 2 = 18.5° and for the 46° halo 
ratio (h 46 ) we use 64=46.5° and 6> 2 =43°. We consider the 
halo features to be present in the phase functions if h 22 
and h 46 are larger than unity and to increase in strength 
with increasing halo ratio. 

3. Results 

Fig. 1 shows the h 22 and h 46 halo ratios as a function of 
crystal aspect ratio and distortion parameter. As expected, 
h 2 2 and h 46 decrease with increasing distortion, and for 
distortion parameter values greater than about 0.15-0.2, 
halo ratios are below unity and thus both halo features are 
not present anymore. Plates show a stronger decrease of 
the h 2 2 ratio with distortion than columns. The h 46 ratio is 
greatest for compact particles and decreases with aspect 
ratio increasingly deviating from unity. Only for crystals 
within a narrow range of distortion parameters between 
about 0.1 and 0.15, the situation occurs that the 22° halo 
feature is present while the 46° halo is generally not, 
which is consistent with previous conclusions by Shcher- 
bakov [19]. 
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Fig. 1. Halo ratios h 2 2 (top) and h 46 (bottom) for hexagonal ice crystals 
with varying aspect ratios and distortion parameters. 
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The general dependence of the halo ratios on aspect 
ratio and distortion parameter is qualitatively consistent 
with the computations shown by Shcherbakov [19]. Since 
Shcherbakov uses an uncorrelated bivariate normal dis- 
tribution [33] to simulate the crystal facet tilt distribution 
[cf. 29] rather than an uniform distribution as used in our 
calculations, the qualitative agreement between our 
results and those of Shcherbakov suggests that the general 
conclusions drawn here are largely independent of the 
assumed crystal facet distribution method. 

Next, we compute h 2 2 and h 46 for distributions of 
distortion levels by applying Gaussian weighting functions 
to the results. For three different aspect ratios, Fig. 2 shows 
the h 2 2 and h 46 halo ratios as a function of varying peak 
and standard deviation values of the Gaussian distribution. 
Note that the weighting functions are set to zero for 
distortion levels outside of the considered range (5=0- 
0.8) and are renormalized. The relative contribution to the 
distributions by relatively pristine crystals (5 <0.15) that 
generally produce halos is given in Fig. 3. 

Fig. 2 shows that the 22° halo feature is very persistent 
within a mixture of pristine and distorted crystals as long 
as the distribution of the distortion parameter is wide 
enough. Perusing Figs. 2 and 3, we conclude that the 22° 
halo is generally present in the scattering phase function if 
the contribution of pristine crystals (5 < 0.15) to the total 
scattering cross section is greater than about 10% in the 
case of compact particles or columns, and greater than 
about 40% for plates. The 46° halo feature, on the other 
hand, is less persistent within a mixture of pristine and 
distorted crystals [cf. 34]. For compact particles, the 46° 
halo feature is only present in the phase function if the 
peak distortion level is below 0.3 and the contribution by 
pristine crystals to the total scattering cross section is 


above about 20%. For plates and columns, the contribution 
by pristine crystals needs to be above about 50% and 70%, 
respectively, before 46° halos appear in the phase function. 

Since rather small contributions (~10%) to the total 
scattering cross section by pristine particles are enough to 
produce 22° halo features, these results potentially explain 
the apparent inconsistency between the frequent sightings 
of 22° halos on one hand and the observed dominance of 
distorted, non-pristine ice crystals on the other hand. 

Given that the 46° halo feature is only present if the 
peak distortion level is below 0.3, while remote sensing 
measurements indicate severely roughened or distorted 
particles with average distortion levels equivalent to 
5 = 0.4 or greater 1 to be prevalent [16,35,17,18,36], our 
results also provide an explanation of the common sight- 
ings of the 22° halo without any detectable 46° halo. This 
explanation can be considered as additional to the expla- 
nation proposed by Shcherbakov [19] based on an aniso- 
tropic roughness of the prismatic facets of hexagonal ice 
crystals that may provide a plausible proxy for ensembles 
of hollow-ended columns and/or bullets. 

Note that the results shown here can be generalized to 
any mixture of particles that produce pronounced halo 
features and those that do not. In the mixtures of ice 
crystals considered here, we assume the crystals that 
produce pronounced halo features to be pristine hexago- 
nal ice prisms with size parameters of the order 100 or 
larger and the particles that do not produce halo features 
to be distorted hexagonal ice crystals in the same size 
range. However, smaller hexagonal crystals and non- 
hexagonal particles also produce no or less pronounced 
22° and 46° halo features [4,5]. Theoretically, the distorted 
crystals in the mixtures assumed in these calculations can 
be replaced by any of such crystals that do not produce 
halo features without altering the general conclusions 
about the contributions of pristine hexagonal crystals with 
size parameters of the order of about 100 or greater 
necessary to produce the halo features. 


4. Summary and conclusions 

Halos at 22° and 46° from the sun are associated with 
pristine, undistorted hexagonal crystals with size para- 
meters of the order of about 100 or greater in ice clouds. 
Although the observability of halos in cirrus clouds may 
depend on cloud parameters such as optical thickness and 
extinction profile, the primary requirement for halos to be 
observable is that the halo features are present in the ice 
crystal scattering phase functions. However, in situ and 
remote sensing measurements indicate the dominance of 
distorted, roughened crystals that produce phase functions 
that generally lack such halo features. Considering that 
plausible mechanisms that could cause distortion and 
surface roughness are stochastic processes such as aggre- 
gation, sublimation, riming, and collisions, ice crystals in 
natural clouds are likely to have considerable distributions 
of varying distortion levels. In this paper, we examined 


1 See Neshyba et al. [27] for the relations between different defini- 
tions of distortion used. 
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Fig. 2. Halo ratios h 2 2 (left) and h 46 (right) for hexagonal ice crystals with distributions of distortion parameters. Colors indicate halo ratio values as a 
function of peak value and standard deviation of the distortion parameter distributions. Results for crystals with aspect ratios of 0.1 (top), 1 (middle) and 10 
(bottom) are shown. 
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Fig. 3. Contribution by relatively pristine crystals (S <0.15) to the total 
scattering cross section as a function of the peak and standard deviation 
of the applied Gaussian distribution function. 


the presence of the 22° and 46° halo features in phase 
functions of mixtures of pristine and distorted hexagonal 
ice crystals. 

We conclude that the 22° halo feature is generally 
present in the scattering phase functions if the contribu- 
tion by pristine crystals with size parameters of the order 
of about 100 or greater to the total scattering cross section 
is greater than about 10% in the case of compact particles 
or columns, and greater than about 40% for plates. For 
compact particles, the 46° halo feature is only present if 
the peak distortion level is below 0.3 and the contribution 
of pristine crystals to the total scattering cross section is 
above 20%. For plates and columns, the contribution by 
pristine crystals need to be above about 50% and 70%, 
respectively, before 46° halo features appear. 

These results potentially explain the apparent incon- 
sistency between the frequent sightings of 22° halos on 
one hand and dominance of distorted, non-pristine ice 
crystals indicated by in situ and remote sensing data on 
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the other hand. Furthermore, the persistence of the 22° 
halo feature in the phase functions and the stronger 
diminishing of the 46° halo with an increasing contribu- 
tion to the total scattering cross section by distorted 
crystals provides a potential explanation of the common 
sighting of the 22° halo without any detectable 46° halo. 
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